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PRKFACE 

The  work  reported  below  is  related  to  forced  vibration  testin'’,  wliicb 
was  performed  on  approximately  one-tbird  scale  model  reinforced  concrete  struc- 
tures. The  tests  were  performed  at  Fort  Polk,  l.ouisiana,  in  the  Summer  and 
Autumn  of  1975  in  connection  with  F.vent  KSSEX-V  by  personnel  of  Waterways 
Experiment  Station,  Vicksburg,  Mississippi.  Those  at  WES  who  have  been  es- 
pecially helpful  in  obtaining  data  needed  to  prepare  this  report  include  Dr.  .1. 
P.  Balsara,  Mr.  Robert  Walker,  Dr.  Sam  Kiger,  Mr.  Roger  Crowson  and  Mr.  Jim 
Ball.  Their  cooperation  is  greatly  appreciated. 

Acknowledgment  is  also  due  Mr.  ,1.  P.  Wright  of  Weidlinger  Associates, 
New  York,  New  York,  whose  suggestions  led  to  pursuing  the  effect  of  backfill 
stiffness  on  radiation  damping.  The  authors  are  also  indebted  to  Professor  H. 

H.  Bleich,  consultant  to  Weidlinger  Associates,  who  pointed  out  the  possibility 
that  friction  at  the  soil-structure  Interface  is  a significant  mechanism  of 
energy  absorption.  The  authors  are  also  grateful  to  Dr.  T.  L.  Geers  of  Lockheed 
Palo  Alto  Research  Laboratory  who  made  available  results  of  analytic  solutions 
for  vibrations  of  a cylindrical  shell  in  an  infinite  elastic  domain. 
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NOTATION 


E , E 
o m 


radius  of  slioll 

trihut.iry  surface  area  associated  with  nodal  point 
wave  speed  (in  >;oneral) 

bar,  shear  and  dilatational  wave  speeds,  respectively 
viscous  friction  coefficient 

depth  of  a point  on  a buried  structure  below  ground  surface 

Young's  m<idulus  of  structure  and  idealized  soil  medium, 
respectively 

matrix  of  friction  forces 
maximum  harmonically-applied  force 

acceleration  due  to  gravity 
thickness  of  structure  wall 
area  moment  of  inertia 

stiffness  of  single  degree-of-f reedom  oscillator 
stiffness  matrix 
width  of  an  element 
mass  matrix 

cycle  number  (subscript) 

typical  nodal  normal  force 

frequency  of  excitation 

matrix  of  externally  applied  forces 

radius  of  arch 


X,  X,  X 


particle  velocity 
maximum  velocity 

weight  density 
weight 

displacement,  velocity,  acceleration 

maximum  displacements  In  successive  cycles  of  oscillation 
maximum  displacement  of  single  degree-of-f reedom  system 
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X 

stat 


6 , y 

s o 

o 

6 

A 
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n 

y 

y . V 

o m 


P 


o ’ 


P 


m 


"h* 


a 

V 


maximum  displacemont  under  statically  applied  load 
non-dimensional  shell  pa;  imeters 

ratio  of  bar  wave  velocity  to  dilatational  wave  velocity 
log  decrement 

angle  subtended  by  two  radii  from  ar  h centerline  to  two  adjacent 
element  centroids 

strain 

fraction  of  critical  damping 
coefficient  of  friction 

Poisson's  ratio  for  the  structure  and  soil,  respectively 
mass  density  of  structure  and  soil,  respectively 
horizontal  and  vertical  components  of  overburden  stress 


O normal  stress  on  surface  of  structure  due  to  overburden 

n 


T maximum  friction  stress 

max 


angle  of  a node  on  the  surface  of  the  arch  (positive  clockwise 
from  horizontal) 

natural  frequency  of  single  degree-of-f reedom  oscillator 


SFXTION  I 
INTRODUCTION 

1.1.  BACKCROUND 

Forced  vibr.it  ion  tests  were  conducted  on  .t  1:3.7  scale  burled  arch 
and  on  two  buried  rectangular  structures  at  the  site  of  ESSKX-V  near  Fort  Polk, 
l.ouisiana.  Tlie  objective  of  the  tests  was  tc  determine  stiffness  and  damping 
cliaracteristics  of  the  structures  before  and  after  backfill  was  placed.  Tliis 
information  is  needed  in  order  to  develop  single  degree-of-f reedom  models  of 
structures  for  use  in  target  evaluation. 

The  geometry  of  the  arch  and  backfill  are  shewn  in  Figure  1-1 . The 

arcli  was  first  tested  in  the  excavation  without  backfill  cover.  Then  b.ickfill 

was  added  and  compacted  (98-100  pcf)  up  to  original  grade  level  and  vibration 

tests  repeated.  All  vibration  testing  on  the  arcti  was  completed  before  t lie  main 

ESSEX-V  event.  The  tests  and  results  are  described  in  detail  in  Ref.  1.  Typical 

results  of  the  tests  are  shown  in  Figure  1-2  in  the  form  of  impedance  (F  /V  ) 

" max  max 

and  quadrature  (F  times  peak  acceleration  times  sine  of  the  phase  angle 
max 

between  them)  versus  frequency.  The  valleys  in  the  impedance  curves  and  the 
peaks  in  the  quadrature  curves  for  the  uncovered  structure  shovr  that  the  peak 
velocity  becomes  large  (resonance)  when  the  driving  frequency  approaches  a 
natural  frequency  of  the  structure.  The  curves  for  the  covered  structure  do 
not  show  such  peaks  or  valleys.  This  is  Interpreted  as  absence  of  resonance. 

The  covered  structure  was  tested  three  times  at  peak  force  levels  of  50,  150 
and  500  pounds,  respectively,  and  the  Impedance  versus  frequency  plots  were 
almost  Identical.  This  result  is  considered  to  confirm  the  absence  of  resonance 
in  the  covered  structure. 

Forced  vibration  tests  were  also  conducted  on  two  1:3.7  scale  rectan- 
gular structures,  which  are  Illustrated  in  Figure  1-3.  Unlike  the  arch,  these 
structures  were  subjected  to  forced  vibration  testing  after  ESSEX-V.  The  orig- 
inal backfill  was  excavated,  the  structures  were  inspected  and  were  found  to  be 
undamaged,  and  sand  backfill  was  placed  by  raining  from  a height  of  5 feet. 
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The  plols  of  impodatno  aiul  <|iiatlraturt‘  vorsus  frequency  for  Lhe  stiffer  struc- 
ture (Structure  30,  l,/li  = A)  with  and  without  backfill,  Figure  1-A,  exhibit 
distinct  resonances.  Similar  plots  for  the  more  flexible  structure  (Struc- 
ture 3B,  h/h  10),  Figure  1-5,  show  resonances  for  the  uncovered  structures. 
Some  evidence  of  weak  resonance  was  also  found  in  the  covered  structure,  but 
this  is  somewhat  a matter  of  interpretation  of  the  data. 

The  observations  that  resonance  did  not  appear  in  the  arch  after 
backfilling  and  that  resonance  was  weak  in  one  rectangular  structure  (I,/h  = 10) 
were  unexpected.  The  results  cast  doubt  on  the  usefulness  of  forced  vibration 
testing  to  develop  parameters  for  SDOF  models  of  buried  structures.  The  results 
also  point  out  a deficiency  in  fundamental  knowledge  of  soil-structure  inter- 
action. 

1.2.  OBJECTIVE 

The  objective  of  the  present  study  is  to  explain  why  resonance  was 
suppressed  in  the  burled  arch  and  to  determine  what  changes,  if  any,  should 
be  made  In  SDOF  models  of  buried  arches.  The  results  of  the  tests  on  the 
buried  rectangular  structure  which  clearly  exhibited  resonance  were  kept  in 
mind  when  various  hypotheses  were  evaluated.  A secondary  objective  was  to 
study  methods  of  estimating  the  equivalent  percentage  of  critical  viscous  damp- 
ing from  the  forced  vibration  test  data  and  also  from  tlie  various  finite  ele- 
ment simulations. 

1.3.  SCOPE 

An  investigation  was  made  of  factors  governing  the  behavior  of  the 
arch  before  and  after  emplacement  of  backfill  using  three-dimensional  elastic 
finite  element  simulation  techniques.  The  scope  of  the  studies  included  ob- 
taining mode  shapes  and  frequencies,  performing  dynamic  modal  analysis  and 
dynamic  step-by-step  analysis. 

The  chief  hypothesis  of  the  analytic  studies  was  that  adding  back- 
fill to  tlie  arch  increased  the  effective  damping  to  the  point  where  resonance 
was  significantly  reduced  or  suppressed  altogether.  This  type  of  effect  occurs 
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in  a single  degree-of-f reedom  oscillator  with  even  modest  amounts  of  viscous 
damping.  Figure  1-h  stiows  that  at  12.3  percent  critical  damping,  forced  vibra- 
tion of  a SDOF  oscillator  at  its  resonant  frequency  produces  only  four  times 
the  di sp l.acement  produced  by  the  same  peak  force  applied  statically;  at  50  per- 
cent critical  damping,  peak  dynamic  and  static  displacements  are  equal  (Kef.  2). 

Possible  sources  of  damping  for  the  burled  arcli  were  identified  as 

fol lows : 

a.  Radiation  damping 

b.  Friction  between  soil  and  structure 

c.  Hysteresis  in  soil  adjacent  to  the  structure 

Structural  damping  in  the  concrete  structure  was  another  possible  source  of 
damping;  it  was  assumed  to  be  small  and  was  therefore  not  considered. 

The  amount  of  radiation  damping  for  each  structure  depends  on  the 
natural  frequencies  of  the  soil-structure  system,  its  geometry  and  the  wave 
speeds  in  the  soil  and  structural  materials.  Differences  in  geometry  between 
the  arch  and  rectangular  structures  and  in  the  respective  backfill  materials 
may  account  for  different  amounts  of  radiation  damping.  This  is  considered 
to  be  a significant  effect. 

Friction  must  be  considered  because  overburden  stresses  (o)  of  2 to 

10  psl  allow  maximum  friction  stresses  (t  ) of  between  .5  psi  to  3 psi  to 

max 

develop,  based  on  an  assumed  coefficient  of  friction  (p)  of  .25  to  .3  (Ref.  3). 
Summing  these  stresses  over  the  surface  area  of  the  structure  indicates  a 
potential  for  energy  absorption  which  is  comparable  to  the  total  of  kinetic 
and  potential  energy  of  the  structure  during  the  forced  vibration  tests.  Since 
the  geometry  of  the  arch  structure  Is  more  conducive  to  generating  motions 
tangent  to  the  soil-structure  interface,  which  enables  the  friction  forces  to 
do  work,  and  since  the  soil  is  probably  in  more  effective  contact  with  the 
arch  than  with  the  vertical  walls  of  the  rectangular  structures,  dissipation 
due  to  friction  may  be  greater  for  the  arch.  Friction  is  therefore  a second 
possible  source  of  damping  which  may  effect  the  arcli  and  rectangular  structures 
to  different  degrees. 


1 


Hysteresis  in  stress-strain  reJations  causes  energy  to  be  absorbed 
by  cycles  of  load-unload  or  reload  which  are  produced  in  the  soil  by  oscilla- 
tion of  the  structure.  The  effective  hysteretic  damping  depends  on  the  amount 
of  energy  absorbed  in  each  cycle.  Although  no  measurements  of  soil  strains  were 
made  during  the  forced  vibration  tests,  the  order  of  magnitude  can  be  estimated 
(Ref.  A)  from 


V 

C = — 
c 


where 


C = representative  strain  parameter  in  the  soil 

V = peak  particle  velocity  in  the  soil  (.05  in. /sec.  to  .5  in. /sec. 
for  500  pounds  maximum  force 

c = wave  speed  in  the  soil  (approximately  12,000  in. /sec.  for  backfill. 

This  leads  to  estimates  of  order  5 10  ^ to  5 x 10  ^ in. /in.  in  a small  volume 

of  soil  adjacent  to  the  structure,  and  smaller  values  at  greater  distances. 

It  is  difficult  to  estimate  the  energy  loss  associated  wltli  these  strain  levels 
because  no  accurate  measurements  were  made  on  the  backfill  material  In  this 
range.  Karthquake-related  research  has  developed  a concept  of  equivalent  crit- 
ical viscous  damping  for  soil  subjected  to  shear  waves  (Refs.  5 and  6).  The 
problem  of  converting  this  information  into  fractions  of  equivalent  critical 
viscous  damping  for  the  present  soil-structure  system  is  formidable  and  beyond 
the  scope  of  the  present  study;  however,  an  upper  bound  can  be  determined  by 
assuming  that  the  structure  experiences  as  much  damping  as  the  adjacent  ele- 
ments of  soil,  which  is  that  associated  with  the  strain  range  given  above. 

This  appears  to  be  about  3 to  4 percent  of  critical  damping. 

The  present  study  concentrates  on  radiation  damping  and  friction. 
Preliminary  estimates  suggested  that  radiation  damping  is  the  dominant  effect 
and  major  emphasis  was  placed  on  investigating  it.  The  final  results  con- 
firmed that  radiation  effects  outweigh  friction  by  4 or  5 to  1. 
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SKCTION  2 

KINITi;  KI.I:MKNT  SIML'I^MION  ok  BKRIKI)  ARCH 

2.1.  INTRODUCTION 

The  first  step  in  simul.it  ing  the  forced  vibr.ition  tests  on  the  .irch 
is  to  obtain  natural  frequencies  and  mode  shapes  associated  witli  deformations 
of  the  structure.  Those  of  the  uncovered  structure  were  compared  with  me.isured 
qu.intities  in  order  to  qualify  tlie  dyn.imic  model.  .Modes  and  frequencies  of  the 
covered  structure  were  compared  with  those  from  tlie  uncovered  structure  to 
indicate  the  extent  of  soil  participation.  .I’ext,  a normal  mode  analysis  was 
performed  using  the  results  of  the  modal  extraction.  The  loading  was  a fre- 
quency sweep  in  whicli  25  cycles  of  a harmonic  load  were  applied  .it  the  irown, 
following  which  the  frequency  of  the  force  was  increased  by  roughly  5 percent 
and  held  for  another  25  cycles.  Undamped  modal  analysis  was  conducted  for 
both  covered  and  uncovered  models,  and  a parametric  study  of  the  effects  of 
modal  damping  on  resonance  in  the  covered  structure  was  made  in  which  12.5, 

25  and  50  percent  critical  damping  in  all  modes  were  assumed.  The  tiiird  step 
was  to  clarify  the  role  of  radiation  damping  by  Investigating  the  properties 
of  tlie  structure  in  a semi- inf  in i te  soil  deposit.  Since  such  a system  does 
not  possess  classical  normal  modes,  the  harmonic  modal  analysis  was  abandoned 
in  favor  of  transient  analysis  using  direct  integration,  and  a non-reflecting 
boundary  condition  of  the  l.ysmer-Kuh lemeyer  type  (Ref.  7)  was  used  at  the  sub- 
surface edges  of  the  soil  Island  in  order  to  simulate  a semi- inf  in i te  soil 
domain.  A parametric  study  was  made  to  study  the  effect  of  variation  in  wave 
speeds  in  soil  on  radiation  damping.  The  final  step  reported  here  is  to  in- 
vestigate the  contribution  of  friction  at  the  soil-structure  interface  to  damp- 
ing of  the  covered  structure.  The  linear  finite  element  code  GENSAP  was  amended 
to  apply  damping  forces  at  interface  nodes.  After  executing  several  check 
problems,  a three-dimensional  transient  analysis  of  the  covered  structure,  in- 
cluding nonlinear  effects  of  friction,  was  performed. 
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2.2.  RKSl'I.TS  OK  MODAL  KXTRACTION 

The  geometry  of  the  finite  element  .simulation  used  throughout  the 
present  study  is  shown  in  Figure  2-1.  A plane  of  symmetry  is  assumed  perpen- 
dicular to  the  long  axis  of  the  structure,  with  the  result  that  only  half  the 
structure  is  represented.  A second  geometric  plane  of  symmetry  exists  which 
lies  along  the  axis  but  this  was  not  assumed  in  tlie  model.  In  order  to  pre- 
serve the  option  to  consider  anti-symmetrical  modes  and  unsymmetrical  loading 
it  was  necessary  to  include  both  sides  of  the  structure.  However,  tliis  option 
was  not  used  during  the  phase  of  the  study  reported  here.  The  finite  element 
discretization  of  the  structure  is  sliown  in  Figure  2-2.  Thin  plate  elements 
(four  node,  linear  curvature  compatible)  are  used  for  tlie  structure.  Hexa- 
hedral  continuum  elements  are  used  for  the  soil. 

The  material  properties  of  the  soil  assumed  in  the  analysis  were  based 
on  secant  moduli  in  uniaxial  strain  from  0 to  1-2  percent  strain  (Kef.  8). 

This  was  based  on  a preliminary  upper  bound  estimate  of  the  stresses  wliich  could 
be  induced  by  the  forced  vibration.  The  properties  derived  in  this  way  are 
shown  in  Table  3-1.  As  is  discussed  below,  subsequent  transient  response  anal- 
ysis showed  that  this  strain  level  is  probably  mucli  too  high.  Altliough  the 
transient  analyses  showed  the  radiation  damping  to  be  sensitive  to  changes  in 
the  assumed  moduli,  previous  numerical  experiments  indicate  that  the  effect  of 
soil  stiffness  on  mode  shapes  and  frequencies  is  much  weaker  (Ref.  9).  A de- 
tailed parameter  study  was  therefore  deemed  necessary,  and  the  modes  reported 
below  are  obtained  on  the  basis  of  lower  bound  soil  moduli. 

The  natural  frequencies  and  mode  shapes  were  obtained  with  a Rayleigh- 
Ritz  procedure  in  which  the  Initial  shape  functions  were  a set  of  harmonic  dis- 
placement functions  suggested  by  the  theory  of  cylindrical  shells.  Modes  of 
the  uncovered  structure  were  extracted  by  placing  the  structure  in  the  excava- 
tion without  any  backfill.  Modes  of  the  covered  structure  were  extracted  by 
adding  soil  elements  having  properties  of  backfill.  Table  2-2  compares  natural 
frequencies  of  the  covered  and  uncovered  models  with  measured  frequencies  of 
the  uncovered  structure.  The  range  of  measured  frequencies  is  somewhat  above 
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tfiose  calciilat  i-d . This  is  partly  because  the  calculations  reveal  lower  modes 
wliich  are  strongly  influenced  by  rigid  body  motion  whereas  the  measurements 
tend  to  miss  tliese.  The  calculations  could  have  been  extended  to  higher  modes, 
but  this  would  have  required  finer  discretization  of  both  structure  and  soil. 
This  was  unnecessary  because  only  a few  def ormat iona 1 modes  are  needed  to 
investigate  the  influence  of  damping  on  respcmse.  The  range  from  mode  5 to 
mode  10  of  the  models  is  rich  in  deformations.  .Measured  mode  shapes  and 
natural  frequencies  for  modes  1 and  2 compared  favorably  with  values  calculated 
for  modes  5 and  7 of  the  finite  element  model  of  the  uncovered  arch.  This  is 
shown  in  Figures  2-2  and  2-3.  Measurements  of  the  deflected  shape  for  mode  3 
were  not  reported  by  WKS.  Hence,  it  is  not  certain  whether  the  mode  shape  is 
the  same  as  that  of  mode  9 of  the  uncovered  analytic  model.  The  calculated 
mode  shapes  of  the  covered  and  uncovered  str\icture  are  given  in  Appendix  I. 

Previous  studies  (Ref.  9)  have  correlated  mode  shapes  of  a structure 
with  and  without  soil  cover  by  summing  the  dot  product  of  eigenvectors  from 
the  two  analyses.  This  correlation  was  part  of  an  effort  to  develop  in  va^uuq 
models  which  would  have  the  same  frequency  characteristics  as  buried  structures. 
The  present  study  investigates  the  damping  associated  with  free  vibration  of 
an  embedded  structure.  This  Information  may  eventually  be  combined  with  the 
previous  data  to  develop  an  ijt  vacuuo  model  with  approximately  the  same  fre- 
quency and  damping  characteristics  ns  the  embedded  structure. 

2.3.  SIMULATION  OF  FORCED  VIBRATION  MODAL  ANALYSIS 

Simulation  of  forced  vibration  at  the  crown  of  both  covered  and  un- 
covered structures  was  performed  using  the  natural  frequencies  and  mode  shapes 
described  above.  The  crown  of  the  arch  was  driven  with  a sinusoidally  varying 
force  whose  peak  magnitude  was  250  pounds  (500  pounds  for  full  structure). 
Starting  at  a frequency  of  75  Hz,  the  frequency  of  the  forcing  function  was 
held  fixed  for  25  cycles.  The  frequency  was  then  Increased  by  5 percent  (to 
78.75  Hz)  and  held  for  another  25  cycles.  This  procedure  was  continued  until 
a frequency  of  141  Hz  was  reached.  Damping  was  assumed  to  be  zero  apart  from 
a small  amount  present  in  the  linear  acceleration  method  used  to  integrate  the 
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mDtial  equations.  The  results  are  shown  in  Figures  2-^  and  2-5  in  terms  of 
velocity-time  histories  at  the  drive  point.  Each  discrete  frequency  and  tiie 
interval  during  wliicli  it  is  applied  is  shown  on  the  time  axis.  Each  natural 
frequency  of  the  finite  element  simulation  is  marked  by  an  arrow  pointing  to 
the  time  interval  (and  hence,  frequency)  where  resonance  miglit  be  expected. 

Figure  2-4  shows  the  ve loci ty-t ime  history  for  the  uncovered  arch 
at  the  crown  (excitation  point).  There  are  three  natural  frequencies  that 
occur  within  the  range  of  the  sweep  and  which  may  appear  as  a resonant  peak. 

The  ones  at  84.2  and  100.5  Hz  (modes  6 and  7)  can  be  recognized  because  of 
the  relatively  high  amplitudes  reached  as  the  forcing  function  sweeps  through 
them.  However,  mode  8 (119.9  Hz)  is  not  detected  at  ail.  This  is  apparentl 
because  it  involves  predominantly  floor  deformations. 

The  corresponding  velocity-time  history  at  the  crown  of  the  covered 
arch,  Figure  2-5,  contains  three  natural  frequencies  within  the  frequency 
range  swept  and  range  of  modes  extracted.  The  first  natural  frequency  in  this 
range  is  83.8  Hz.  However,  this  is  the  same  floor  mode  that  did  not  exhibit 
resonance  at  the  crown  of  the  uncovered  arch  and  it  is  believed  that  it  goes 
undetected  again  for  the  same  reasons.  However,  modes  9 and  10  are  observed. 
These  are  at  100.5  and  110.8  Hz.  Because  of  their  close  spacing  extra  care 
must  be  taken  in  examining  and  Interpreting  the  results. 

The  occurrence  of  resonance  in  both  models  indicates  that  much  greater 
damping  was  needed  to  simulate  the  physical  experiment.  In  order  to  estimate 
the  equivalent  viscous  damping  required  to  suppress  resonance,  which  would  agree 
qualitatively  with  the  observations  on  the  covered  arch,  three  modal  analyses 
were  made  In  which  the  percentage  of  critical  damping  was  12.5,  25  and  50  per- 
cent. In  each  case  the  damping  was  assumed  to  be  constant  in  all  modes.  The 
calculations  are  thus  similar  to  the  one  whose  results  are  shown  in  Figure  2-4, 
except  that  the  damping  is  greater.  The  results  sliow  that  resonance  occurs 
at  12.5  percent  and  that  it  Is  suppressed  at  50  percent.  At  25  percent  critical 
viscous  damping  there  are  only  slight  remnants  of  resonance.  There  is  motion 
in  both  the  25  and  50  percent  cases,  but  the  amplitude  is,  to  first  order,  inde- 
pendent of  frequency  of  the  forcing  function. 
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Ttif  result  described  above  indicates  tli.it  .in  eltective  damping  of 
20  to  10  percent  is  required  to  bring  the  an.ilytic  model  i.Uo  qualitative  agree- 
ment with  the  test  data.  The  questions  to  be  .iddres-ied  next  are:  Wli.it  are  the 
sources  of  damping  and  is  their  aggregate  etlect  gre.it  enougli  to  suppress  reso- 
nance? 

2.4.  TRj'\NSIKNT  ANALYSIS  OK  COVKRHO  AKCH 

Three  damping  mechanisms  have  been  postul.ited  in  Sc-ction  1 as  con- 
tributing to  the  total  d.amping  of  tlie  arch: 

a.  Radiation  damping 

b.  Friction 

c.  Hysteresis  in  backfill  material 

Radiation  damping  is  associated  with  wave  propagation  in  an  infinite 
or  semi-infinite  domain  and  is  not  correctly  represented  in  a bounded  domain 
such  as  the  one  analyzed  by  modal  extraction  methods.  To  represent  a semi- 
infinite soil  domain,  a Lysmer-Kuhlemeyer  non-reflecting  boundary  condition 
(Ref.  7)  was  apjlied  to  the  subsurface  boundaries  of  the  soil  Island.  This 
model  no  longer  possesses  classical  normal  modes,  and  hence  resort  is  made 
to  a direct,  step-by-step  integration  method.  The  possibility  that  friction 
between  the  soil  and  structure  may  absorb  a significant  amount  of  vibrational 
energy  arises  from  the  significant  confining  stresses  (2-10  psi,  depending 
on  depth)  due  to  overburden.  Sliding  friction  introduces  nonlinearity  into 
the  system  of  equations,  which  is  another  reason  for  using  direct,  step-by-step 
Integration  methods.  As  is  explained  above,  soil  hysteresis  is  omitted  for 
lack  of  experimental  data  on  which  to  base  a constitutive  model  and  structural 
damping  is  omitted  because  it  is  small. 

In  the  following  analysis,  damping  is  investigated  by  observing  decay 
in  the  amplitude  of  free  vibration.  Motion  is  induced  by  applying  a triangular 
pulse  whose  duration  is  about  equal  to  .005. second  (approximately  the  half 
period  of  the  100.7  Hz  mode).  The  total  duration  of  the  calculation  is  .108 
second,  or  about  21  times  the  duration  of  the  load;  case  2 is  an  exception  due 
to  being  terminated  prematurely  by  the  computer  center  operator. 
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2.4.1.  KKKKCT  OF  SOIL  STIKFNKSS  ON  10\1)IAT10N  DA.MPING 

Ciilciilat  ions  were  m.ule  to  Invest  iK-itt  radiation  damping  and  the  in- 
fluence of  wave  speed  in  the  backfill  and  surrounding  soil  on  the  amount  of 
radiation  damping.  The  finite  element  simulation  showed  that  the  amount  of 
radiation  damping  is  sensitive  to  variations  in  the  assumed  wave  speeds  wliich 
are  within  the  range  of  experimental  uncertainty.  All  of  the  analyses  were 
based  on  measurements  of  uniaxial  stress-strain  properties  performed  on  undis- 
turbed samples  taken  from  the  LSSEX  test  site.  The  site  profile  and  repre- 
sentative stress-strain  data  are  given  in  Ref.  8.  Among  the  most  Important 
data  is  the  uniaxial  stress-strain  curve  for  the  tamped  backfill  surrounding 
the  arch,  which  is  shown  in  Figure  2-6. 

In  Interpreting  the  data  for  use  in  the  calculations,  attention  was 
given  to  possible  differences  between  the  actual  stress-strain  curves  for 
backfill  in  the  field  and  the  representative  curve  given  in  Figure  2-6.  There 
is  the  possibility  of  disturbance  during  sampling;  the  effects  of  sampling 
procedures  on  moduli  are  not  well  known  at  the  exceptionally  low  stress  levels 
of  interest  in  the  present  study.  This  problem  is  reflected  in  the  observa- 
tion that  ^ situ  seismic  wave  speeds  are  In  some  cases  (Ref.  10)  hlglier  tlian 
laboratory  stress-strain  measurements  would  lead  one  to  expect.  Also,  there 
is  variability  in  such  properties  as  the  extent  of  compaction  of  the  backfill, 
wiiich  is  reflected  neither  in  the  stress-strain  curve  in  Figure  3-5  nor  in  the 
calculat Ion. 

A second  question  is  whether  the  secant  or  the  initial  tangent  moduli 
of  the  experimental  stress-strain  curves  should  be  used  as  a basis  for  the 
elastic  moduli  of  the  finite  element  model.  A parametric  <tudy  was  performed 
to  investigate  how  backfill  and  in  situ  soil  stiffness  influence  'adlation 
damping.  A lower  bound  on  moduli  in  tiie  backfill  and  other  soil  layers  was 
established  on  the  basis  of  the  secant  modulus  from  0 to  1-2  percent  strain. 

Tlie  secant  modulus  in  the  backfill  is  about  6,700  psi,  which  leads  to  a dilita- 
tional  wave  speed  of  about  533  psi.  The  wave  speeds  assumed  for  the  backfill 
and  the  other  soil  layers  are  summarized  in  Table  2-3.  The  results  of  the 
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I'^iKiilat  ion  using  these  properties  are  illustrated  by  a ve  loc  i ty-t  ime  history 
at  tiie  crown  in  Figure  2-7.  The  response  is  lightly  damped  and  Is  unlike  tlie 
response  of  tlie  physical  arcli.  Suspecting  that  tliis  is  due  to  ineffective 
transfer  of  energy  in  deformat lonal  modes  to  tlie  soil,  a second  calculation 
was  made  in  wliich  the  soil  properties  were  based  on  the  Initial  tangent 
moduli,  which  are  about  four  times  stiffer.  The  justification  for  this  is 
that  the  stress  levels  developed  in  the  backfill  during  the  tests  were  prob- 
ably of  order  0.1  psl  (order  .01  bars),  so  that  the  initial  tangent  modulus 
is  more  representative  of  the  stiffness  governing  radiation  damping  effects 
in  the  forced  vibration  tests  (thougli  not  necessarily  under  weapons  effects 
loading)  than  is  the  secant  modulus.  The  results  of  the  second  calculation, 
shown  as  ve loc i ty-t ime  history  of  the  crown  in  Figure  2-8,  indicate  a sub- 
stantial increase  in  radiation  damping.  Since  this  result  is  central  to 
the  p'^-pose  of  the  present  study,  it  was  decided  to  perform  another  calcula- 
tion which  winild  establish  a reasonable  upper  limit  on  the  radiation  damping. 

In  this  calculation,  the  initial  tangent  moduli  were  doubled.  This  is  a 
modest  increase  in  the  recommended  values  which  in  the  opinion  of  the  authors 
could  have  occurred.  The  results,  shown  in  Figure  2-9,  confirm  the  trend 
toward  increasing  radiation  damping  with  Increasing  wave  speeds  in  the  soil. 

It  appears  that,  as  the  impedance  of  the  soil  approaclies  tliat  of  the 
concrete,  coupling  Improves  between  deformational  modes  of  the  structure  and 
the  soil.  In  the  extreme  case  of  a perfect  Impedance  match  between  structure 
and  medium,  the  radiation  damping  would  be  very  high;  in  the  opposite  extreme 
case  of  a structure  j_n  vacuuo , there  would  be  no  damping.  In  case  3 (2  X tan- 
gent moduli)  where  the  Impedance  ratio  of  concrete  to  backfill  is  about  1 to 
15,  the  damping  is  estimated  to  be  the  equivalent  of  12  to  21  percent  critical 
viscous  damping,  depending  on  Interpretation  of  tlie  response  of  the  crown. 

In  order  to  confirm  the  trend  obtained  from  the  finite  element  calculations, 

Dr.  T.  L,  Geers  of  Lockheed  Palo  Alto  Research  I.aboratory  was  asked  to  compute 
analytically  the  equivalent  fraction  of  critical  viscous  damping  for  an  Infinite 
cylinder  having  roughly  the  same  geometric  and  material  properties  as  the  arch 
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and  which  is  embedded  in  an  infinite  medium  liavlng  roughly  the  same  properties 
as  those  assumed  for  the  backfill  in  the  finite  element  calculations.  fn  spite 
of  significant  differences  between  the  arch  and  the  cylinder,  the  trends  sub- 
stantially agree.  Results  of  the  cylinder  analysis  are  shown  in  Appendix  11. 

2.4.2.  EFFECT  OF  FKICTION  ON  DAMP  I NO 

The  energy  of  vibration  which  is  absorbed  by  friction  between  the 
soil  and  the  surface  of  the  structure  was  investigated.  An  ideal  computer 
simulation  of  friction  would  include  three-dimensional  geometry,  non-reflecting 
boundaries  and  the  capabiiity  for  slip  at  the  soil-structure  Interface  and  a 
Coulomb-type  friction  slip  criterion  in  which  the  direction  and  magnitude  of 
the  friction  force  depend  on  the  relative  velocity  between  the  soil  and  struc- 
ture. A step  was  taken  toward  this  ideal  simulation  by  modifying  the  existing 
GENSAP  model  to  Include  frictional  forces  tangential  to  the  interface  and  in 
the  direction  opposite  to  the  absolute  velocity  of  the  structure.  It  is  be- 
lieved that  this  simulation  is  adequate  for  purposes  of  this  study.  The  maxi- 
mum frictional  stress  which  can  be  applied  was  assumed  to  be  proportional  to 
the  normal  stress  acting  on  the  structure,  which  depends  on  the  depth  of  each 
nodal  point  below  the  surface.  The  relationship  between  the  overburden  and 
frictional  stresses  is  shown  in  Figure  2-10.  The  frictional  shear  stress  was 
multiplied  by  the  tributary  area  for  each  node  and  the  resulting  force  was  in- 
cluded in  the  discretized  equation  of  motion  as  follows: 

Mx  + Kx  - P - F 


where 


friction 
ever  the 


M,  K “ global  mass  and  stiffness  matrices 
X,  X “ acceleration,  displacement  vectors,  respectively 
P “ externally  applied  load  (forcing  function) 

F “ friction  force  (applied  only  to  soil-structure  interface  nodes) 

An  initial  attempt  was  made  to  compute  F according  to  a simple  Coulomb 
law  In  which  the  maximum  magnitude  of  friction  force  was  applied  when- 
absolute  velocity  of  a node  on  the  structure  exceeded  a threshold  value 
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The  direction  of  the  force  was  assumed  always  to  be  opposite  to  the  velocity 
of  the  structure.  Tliis  attempt  was  unsuccessful  because  the  system  was  too 
strongly  nonlinear  for  the  integration  algorithm  to  handle  it  correctly.  The 
results  depended  on  the  time  step  used.  Consequently  a nonlinear  viscous 
friction  model,  in  which  the  frictional  force  was  assumed  to  be  proportional 
to  velocity,  was  used. 

The  nonlinear  viscous  friction  model  used  was  adapted  from  the  work 
of  Den  Hartog  (Ref.  2);  a viscous  single  degree-of-f reedom  model  in  which  the 
work  per  cycle  done  by  the  friction  force  is  equivalent  to  a constant  force 
Coulomb  model.  In  the  viscous  model,  it  is  assumed  that 

F = Cx  (2-1) 

where 


and 


triox 


o 


(2-2) 


f » peak  friction  force 
w * frequency  of  the  SDOF  system 

= peak  displacement  of  the  SDOF  system 

This  model  exaggerates  the  effects  of  friction  by  assuming  that  die  f r i ct ion  forces  are 
related  to  the  absolute  rather  than  to  the  relative  velocity.  It  is  believed 
that,  because  the  loading  Is  from  Inside  the  structure,  the  absolute  motion  of 
the  structure  is  reasonably  close  to  the  relative  motion.  This  would  not  neces- 
sarily be  the  case  for  engulfment  by  a shock  wave. 

The  implementation  of  this  type  of  model  in  GENSAP  required  numerous 
coding  changes,  and  two  check  problems  were  devised  to  test  the  new  model. 

These  are  described  in  Appendix  111. 

An  analysis  of  the  covered  arch  was  performed  which  was  aimed  at  pro- 
ducing an  upper  bound  estimate  for  damping  Including  the  effects  of  friction 
and  radiation  damping.  The  friction  parameter  f In  Equation  2-2  was  calculated 
for  each  node  as  follows: 
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f = u o A (2-3) 

n 

where 

a - iiuMual  stress  at  each  surface  node  due  to  overburden 
n 

A = tributary  surface  area  associated  with  the  node 
U = coefficient  of  friction 

The  nodal  friction  force  is  applied  tangentially  to  the  arch  in  the  direction 
opposite  to  the  tangential  component  of  velocity.  It  is  the  large  cumulative 
area  which  magnifies  the  friction  force  and  wliich  makes  friction  potentially  a 
significant  effect.  In  the  present  example  p = 0.4  was  assumed  on  the  basis 
of  data  presented  in  Ref.  3. 

The  wave  speeds  in  the  soil  are  based  on  two  times  the  initial  tan- 
gent moduli  (see  Table  2-3);  hence,  radiation  damping  is  maximum  for  the  range 
of  parameters  considered  in  this  study.  This  is  not  a situation  which  maximizes 
dissipation  due  to  friction  because  large  radiation  damping  reduces  velocities 
which  in  turn  reduce  friction  forces.  However,  the  total  damping  is  maximum, 
or  nearly  so. 

The  results  of  tlie  calculation  are  best  Interpreted  when  compared 
with  the  response  without  friction,  which  Is  done  in  Figure  2-li  for  the  crown. 
Friction  reduces  the  amplitudes  of  oscillations  at  intermediate  and  late  times, 
but  the  effect  is  not  dramatic.  The  response  of  the  model  also  contains  higher 
frequency  motions  at  intermediate  and  late  times.  This  Is  due  to  the  assump- 
tion that  the  friction  force  is  proportional  to  the  absolute  velocity  and  has 
the  opposite  direction.  High  frequency  oscillations  are  even  more  characteristic 
of  Coulomb  than  of  viscous  damping  (Ref.  2).  Another  aspect  of  the  response 
with  friction  is  that  the  first  downward  motion  Is  greater  than  in  the  case 
without  friction.  The  reason  for  this  is  apparently  discretization  error  in 
space,  time  or  both.  The  frictional  forces  are  applied  tangentially  to  the 
true  curve  of  the  arch,  which  is  only  approximated  at  each  node  by  the  two  inter- 
secting flat  plate  elements.  Thus,  the  lines  of  action  of  the  internal  resist- 
ing forces  are  not  exactly  correct,  and  this  error  is  compensated  in  the  inertia 
terms  In  order  to  achieve  equilibrium.  Also,  at  early  times  the  friction  forces 
are  the  same  order  as  the  inertia  forces,  which  means  that  nonlinearity  is 
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strong.  The  integration  algorithm  is  not  capable  of  handling  strongly  non- 
linear systems  perfectly,  and  some  errors  undoubtedly  occur.  For  the  main  part 
of  the  calculation,  the  friction  forces  are  smaller  than  the  Inertia  forces  and 
the  conditions  for  accurate  integration  are  more  nearly  met. 

2.5.  EVALUATION  OF  PERCENTAGE  OF  EQUIVALENT  CRITICAL  VISCOUS  DAMPING  FROM 
CALCUIJVTIONS 

Evaluation  of  an  equivalent  fraction  of  critical  viscous  damping  from 
the  results  of  free  vibration  computations  is  necessary  to  decide  whether  the 
damping  mechanisms  considered  analytically  can  account  for  the  suppression  of 
resonance  in  the  physical  vibration  te.sts.  It  is  difficult  to  arrive  at  a 
clear  and  unambiguous  fraction  of  critical  damping  because  the  response  of  the 
arch  Involves  several  modes.  It  is  therefore  hard  to  apply  the  conventional 
log  decrement  Interpretation  to  the  di-splacement-t ime  history  because  it  is 
difficult  to  Identify  successive  maxima  in  the  same  mode.  It  is  also  probable 
that  the  fraction  of  critical  damping  varies  from  mode  to  mode,  as  is  shown  by 
the  results  given  in  Appendix  II.  In  the  following,  the  log  decrement  approach 
Is  applied  to  successive  maxima  according  to  the  expression 

-Hn  - 6 ' — (2-0 

^n 

where 

X , X . , « successive  maximum  amplitudes 
n n + 1 

r|  “ fraction  of  critical  viscous  damping 
6 “ log  decrement 

Pairs  of  successive  maxima  are  obtained  from  each  calculation  and  Equiitlon  2-4 
is  solved  for  n-  To  try  to  reduce  bias  in  selecting  maxima  for  each  case,  two 
pairs  of  maxima  are  chosen  wherever  possible.  The  pairs  selected  here  are  the 
first  and  second  and  the  second  and  third. 

Displacement-time  histories  at  the  crown  from  the  four  cases  considered 
are  shown  in  Figures  2-12a  through  2-12d.  In  each  case,  the  points  chosen  for 
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successive  maxima  arc  indicated  by  the  numbers  1,  2 and,  in  must  cases,  '3. 
Table  2-4  shows  the  fractions  of  critical  damping  for  each  case.  The  trend  is 
in  the  direction  of  greater  damping  for  stiffer  soils;  it  is  increased  further 
by  accounting  for  friction.  However,  not  all  damping  values  are  consistent 
with  this  pattern;  it  is  necessary  to  select  the  values  marked  by  an  asterisk 
to  obtain  tlie  trend  described  above.  However,  visual  inspection  of  the  dis- 
placement-time histories  supports  the  interpretation  that  the  cases  are  ranked 
in  Table  2-4  in  order  of  decreasing  damping. 
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Figure  2-1.  Discretization  for  Finite  Element  Mesh. 
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Figure  2-12a.  Displacement-Tine  History  at  Crov.'n  of  Buried  Arch 
on  Secant  Moduli. 
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TAUl.K  2-  5. 

SOU.  I'KOI’I'KTIKS  I'SHD  IN  I’ARA.NK'I'RI  ('  STUDY  ()[■'  RADIATION 
AKSOCIATKD  WITH  KRKF  VIlIRyVIION  OF  RL'RIKD  ARCH 


Dfjitli  (Ft.)  (fXT“) 


SiH'.itU  Miiduli 

( 1-2 ■ Sf  rain  ) 

l.ayiT  1 
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') 
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Bark  fill 
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lan^imt  Moduli 

l.ayrr  1 
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TAlll.K  2-4. 

I'KACTIONS  OK  CKITICAI,  VKSCOUS  DAMP  I NO  KKOM 

KRKl.  VIHRATION  ANAl.YSKS  OK  Rl'RIKlJ  ARCH 

Case  ^“2  Oyclos  2-3  Avera^f 

Soi  1 pro|H  rt  ii's  based  on 
2 X tan>;iMU  moduli;  fric- 


t ion  included 

. lA 

.26 

20 

.Sc.il  properties  h.ised  on 

2 X tangent  moduli 

. 12 

.21* 

17 

Soil  proj)erties 
tangent  modul  i 

b.ised  on 

. 1 7* 

.Oh 

12 

Soil  properties 
Secant  moduli 

based  on 

.12* 

*If  tiK-sf  case.s  are  .si-lected  ns  bcinp,  most  ri‘()rfsonI  at  i ve , a trend  ol  increasing 
damplnK  with  increasing  soil  stiffness  is  obtaiiud.  including  friction  increases 
d.imping  to  the  level  which  suppressed  resonance  in  the  1 requency  swi'cp  simulation 
using  modal  analysis. 
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SKCTION  5 

INTr.Rl’KKTATION  OF  FXPFR  IMFNTAl.  DATA  AXD  ( ORRF.I.AT  rON  WIIM  ANALYSIS 

1.1.  INTRODK  TION 

The  piirpi'se  of  this  sei  t ion  is  to  discuss  the  expo  r i n-i  lit  .1 1 nic.isiin-- 
nionts  on  the  ;irc)i  in  light  of  the  nn.llvsos  described  in  the  privioiis  section 
• md  to  explain  wliy  at  least  one  of  tlie  buried  rc.“C  t angii  1 ar  striictnres  exiiihited 
strong  rc-sonance  whereas  in  tlie  buried  arch,  resonance  was  suppressed. 

1.2.  KSTI.t'WTK  or  DAMl’lNt;  iN  TKSTS  ON  HL'RIKI)  ARCH  AND  tiOMi’AR  1 SON  Wnil  ( Al.t  I LATIONS 

One  mi  t hod  ot  e tiniating  the  e<iuivalent  traction  c'f  criti  -.il  viscous 
damping  of  the  covered  .irch  from  experimental  data  is  to  regard  it  as  a single 
degree-c'f- ri-edicm  I’Srillator,  Although  thisv'is  ni^t  the  case  it  is  a reasonable 
idea  1 i z;it  ion  because  the  I ri-cjuency  sweep  tends  to  excite  each  normal  nii'de  in 
turn  and  that  mode  temporarily  dominates  ail  otliers.  lor  st  isuiy  state  forced 
vibration  i>t  a single  di'gne-o  f - 1 rc-edom  oscillati'r  we  have; 

/L 

- (1-1) 
p / . • . 

max  > ( 1 - p'  / . ■ )’  + 4'  p' 


where 

V = peak  vi'locitv 

max 

F = peak  harmonic  exciting  force 

max 

p = frecpiency  of  excitation 
la  = natural  t recpic-ncy  of  system 
k = static  stiffness 
1,  * fraction  of  critical  damping 

If  the  arch  is  siihiected  to  exiitation  of  its  natural  frecpiency  (tor  example, 
p * , ' IK)  Hz),  the  impedance,  taken  from  Figure  1-2,  is  1.5  ‘ in'*  1 h . -sec . / in . 
(vibrations  in-phase)  or  ').()  » U)''  lb. -sec. /in.  (vibrations  IHo'^  out -of -phase) . 
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rlif  St. It  it-  .'lillni'ss  i>!  I lie  .irili,  defim-d  hy  the  furee  reijiiired  te 
prediit'e  ;i  unit  iletUnlion  at  the  irnwii,  is  k ' 5 ■ If)^  Ih./tl. 

Simplitvin>;  hqn.ition  4-1  we  Dht.iin 

2K  . ,/V 

fTiei  A rnd X / •»  > » 

’'  = ■ k • >' 

Siihst  i t lit  ini;  exner  imeii  t a 1 1 v measured  values  ut'  I /V  inte  this  exp  re.ss  i on , 
the  ri'sults  presented  in  Tahle  J- 1 are  obtained.  K(|uili()n  i-2  has  hein  evaliiati'd 
• It  p = I(l')  llz  ami  222  Hz  because  resonance  is  observed  in  tlie  uncovered  .ircii  in 
this  ran^;i'  and  because  p = 105  Hz  is  in  the  center  of  the  ran^e  considered  by 
the  mod.il  analysis.  Section  2.').  wliicli  indicated  that  25  percent  critical  damp- 
ing is  tlie  minimum  needed  to  suppress  resonance.  I'he  results  presenteil  in  the 
table  are  not  reason.ible  except  in  t be  case  of  vibr.iLors  in  phase  at  i'  = 105  Hz 
when-  the  uncovered  structure  with  only  6 percent  critical  d.iniping  would  exhibit 
resonance  whereas  the  covered  structure  with  41  percent  would  not.  The  other 
values  are  apparently  too  high.  This  exercise  illustrates  the  ilitficulty  of 
evaluating  the  fractions  of  critical  damping  from  the  experiments.  It  is  only 
slightly  easier  to  estimate  the  cflective  damping  in  the  calculations,  as  is 
desiribed  in  Section  2.5. 

Ihe  failure  of  this  .ipproach  to  yield  reasonable  values  for  damping 
in  the  experiments  was  unfortunate  since  it  metint  that  there  was  no  method  of 
comparing  quan  t i t at  i vi’ I y the  measured  and  calculated  tractions  ot  iritical  liamp- 
ing.  Instead,  it  must  be  argued  that  .since  the  modal  analy.sis  indicate.s  that 
about  25  percent  critic.il  damping  suppresses  resonance  .ind  since  this  amount  can 
be  ai  counted  tor  in  the  finite  element  simulation  it  w.ive  siieeds  are  based  on 
2 X tangent  moduli,  it  is  possible  that  the  fundamental  phenomena  observed  in 
the  tests  .ire  accounted  lor  in  the  calculation. 

1.3.  1)1  l•■H■;RI■:NT  RK.SPONSKS  l)K  ARCH  A.M)  RKCTANClil.AR  STRUCTHRKS 

The  buried  rect  .angu  1 ;ir  structure  II)  (l./b  = 4)  clearly  e.xhibited 
resonance  in  thi'  experiments  .is  shown  in  riguri'  1-4.  The  buried  rectangular 
'.t  riicturi-  lb  (l./h  « 10)  less  clearly  exhiliited  resonance  in  the  experiments 
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IS  shown  in  KiRuro  l-h.  but  WKS  i-onsidiTs  llint  t licro  w.is  wt  .ik  rosonnin  o ovor  t 
bio.ul,  low  t fi'qm-iKV  r.inRf,  Rof.  1.  In  contrast,  the  arch  structure  is  con- 
siiieri'd  to  have  exhibiti-d  no  ri-sonance.  A major  factor  i-on  t r i but  i hr  to  hiRher 
d.impiiiR  for  the  arch  t h.in  lor  the  rectaiiRular  structures  is  dilferent  baikfill 
I'onditions.  Thi'  .arch  structure  was  placed  in  native  backfill  which  was  com- 
pacted civnamical  ly  in  stages.  The  resulting  material  had  a st  ress-st  r.iin  curve 
which  was  initially  concave  to  the  strain  axis  and  which,  as  a result,  h.id  a 
high  seismic  wave  veloc  ity.  The  calculations  described  above  have  clearly  demon- 
strated that  this  situatic)!!  leads  to  strong  coupling  between  structure  and  the 
backfill  and  surrounding  soil  and  to  relatively  high  radiation  damping.  The- 
rectangular  structures  were  excavated  alter  the  m;i  i n T.SSKX-V  event,  and  new  s;md 
backfill  was  rained  into  place.  The  properties  of  the  sand  backfill  were  not 
measured  directly,  but  they  are  considered  by  WKS  to  be  about  the  same  .as  Cook's 
H.ivou  s.md.  fhi'.  sand  has  a stress-strain  curve  which  is  convex  to  the'  strain 
■ ixis  and  has,  as  .1  result,  .1  relatively  low  seismic  wave  velocity.  The  uniaxial 
st  ress-st  r.i  i n curve  for  Cook's  H.iyou  sand  is  shown  in  figure  2-b  along  with  the 
initi.il  tangent  modulus  (measured)  anci  2 X in  rial  t.angent  modulus  for  the  arch 
backfill.  The  initial  t.ingent  modulus  for  Cook's  Bayou  sand  is  .ibout  the  same 
as  the  sec  , ant  modulus  used  in  c-ase  1 tor  the  arch.  Table  2-1,  which  produced 
small  values  of  radiation  damping.  It  is  also  possible  th.it  the  procedure  of 
raining  b.ickf i 1 1 alongside  .1  structure  with  vertical  walls,  such  as  the  rectan- 
gul.ir  structures,  mav  not  produce  strong  coupling  between  b.ackf  i 1 1 and  structure 
during  vibration  testing. 

To  summarize  tin-  different  backfill  conditions,  the  arch  was  well 
j coupli'd  to  a rc'latively  high  wave  speed  backfill  where. is  the  rect.mgiil.ir  st  ruc- 

J tores  wc're  emplaced  in  lower  wave  speed  backfill.  The  net  result  was  that  radia- 

i t ion  damping  was  certainly  higher  in  the  arch  than  in  the  rect.ingiilar  structures, 

I and  possibly  it  was  much  higher. 

I 
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TAlil.H  j-1. 

rRi\c:TIONS  OK  critical  D/VMIMNC  LSTIM/VIKI)  I-'ROM  'IIIK 
ICXI’FRIMKNTAI.  DATA  ON  ARCH  USING  KQUATION  3-1 


Ci'vored 

Uncovered 

In  Phase 

Oul  of  Phase 

In  Phase 

Out  of  Phase 

.41 

.81 

.06 

.27 

2.2 

4.2 

.55 

1.  I 
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siN(;i.K  I)[:(;ri.i;-oi'-khi:i:i)om  taiwmi  ti-ks  i-kom  forci.d  vi  krai  ion  ii  sts 

Tlif  value  ol  forcail  vibratiun  tests  in  deve  lup  i ng  parameters  tor 
b i >;h  1 y simplified  structural  models,  stub  as  .SlK)!'  oscillators,  aitpear.  limited 
in  liRiit  of  experimental  and  analytic  fintiiiiK^-  These  tests  are  probabl-.-  .i 
Kood  w.iy  to  obtain  data  on  frequencies  in  the  elastic  reK*™’  1<'"’K  as  tlie 
fretpiencies  show  up  (unlike  tiie  arch  case)  and  so  ionn  as  I lu'  backfill  of  tlie 
model  is  tamped  simil.iriy  to  that  of  the  prototype.  file  |)ossibility  of  les-.ened 
participation  ot  soil  mass  due  to  loose  backfill  must  be  recottnixed  Itotb  for 
test  and  target  or  prottjtype  it  met  tires.  fable  2-2  shows  that  the  fretpu-ncies 
of  covered  and  tincovered  stru  tures  differ  significantly;  tliese  represent 
lower  .tnd  upper  liounds,  respectively,  on  the  el  feet  of  soil  mass  partiiipa- 
tion.  Ti’.e  analvst  who  intends  to  use  forced  vibration  test  dat.i  mtc.t  eithir 
know  where  the  backfill  conditions  of  his  target  structtire  I'lace  it  in  the 
freqtiency  range  t.hown  in  fable  2-2,  or  be  mtist  be  prepared  to  asstime  upper 
.inti  lower  boiintis  tor  the  f retpiency  of  his  motiel  .iiul  contluct  at  least  two  anal- 
yses in  order  to  bound  tin-  response. 

■fbe  test  data  .and  analysis  a i so  show  th.at  the  contribution  ot  ridi.i- 
tion  to  tlie  total  d.imping  tlepetuis  s i gn  i f i cant  1 y on  the  wave  speeds  in  the  soil. 
Here  again  the  backfill  stillness  in  the  test  must  bi  like  that  of  the  proto- 
type to  be  sure  that  this  I'ffect  is  properly  .iccountid  for.  Since  forced  vibra- 
tion tests  generate  stresses  much  lower  tlian  w»Mpons  efiects  lo.iding,  tin  effec- 
tive wave  speeds  in  the  test  are  eitlu-r  higher  or  lower  t h.in  in  the  prototype 
depending  on  whether  the  uni.ixial  stress-strain  curve  is  st it  fining  or  sol  ton- 
ing in  the  range  ot  interest. 

Anottier  point  related  to  different  stress  levels  in  vibration  tests 
and  weapons  effects  loading  is  the  effect  of  hysteresis  in  soil  s t ri'ss-s t ra i n 
relations.  The  role  of  liysteresis  in  damping  structural  vilirations  is  .at  pre- 
sent unknown  and  should  bi'  investigated.  If  hysteresis  were  s i gn i I leant,  it 
would  be  another  deficiency  of  the  forced  vibration  ti'sts. 


Oni'  w.iv  in  wliicli  ti'sts  like  tlinsc  <'n  t Ik-  arcli  .iml  t)nx  s t rin  t tirf; 
i-niild  bo  niiidi'  ninri'  valuablo  is  to  oonoontrato  mi'ro  lioavily  l lian  is  done 
ctirn-iuly  on  transionl,  tree  vibration  bobavior.  Tin's  would  rivo  anotiur, 
porliaps  rloaror,  pioturo  of  dampinj;. 
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SKCTION  '■> 

Sl'M>WKY  AND  rONCLUS  IONS 

A finite  element  s imii  1 .it  ion  of  the  forced  vihr.ition  testa  on  the 
biiriicl  .ircti  strucliiro  ,it  I',SSKX-V  was  dove  loped . Tlie  model  of  Lite  structure 
surrounded  by  a iinite  amount  of  soil  was  qualified  by  calculating  natural 
frequencies  and  mode  shapes  and,  for  the  uncovered  structure,  comparing  them 
with  measurements.  Modal  analysis  of  the  so i 1 -st rue t ure  system  showed  that 
drive  point  motions  increased  dramatical Iv  in  both  covered  and  uncovered  cases 
when  the  I requency  of  the  oscillating  force  at  the  crown  approached  a natural 
fri'quency  of  the  system.  Modal  analysis  of  the  covered  structure  was  repeated 
using  different  assumed  fractions  of  critical  viscous  damping.  The  resonant 
behavior  of  the  structure  was  suppressed  when  about  2'j  percent  critical  damp- 
ing in  all  modes  was  assumed. 

Railiation  damping,  friction,  hysteresis  in  the  surrounding  nuulium 
.and  stivatural  damping  were  identified  as  possible  mechanisms  for  absorbing 
the  energy  of  vibration  in  the  physical  test.  To  evaluate  the  contributions 
of  r.idiation  damjiing  and  friction,  transient,  free  vibration  of  the  buried 
structure  was  simulated  using  a direct  integration  approach.  Nonre f 1 ec t ing 
boundary  conditions  were  used  on  the  subsurface  boundaries  of  the  soil  island 
to  simulate  a semi- inf  in i te  soil  domain.  It  was  found  that  the  amount  ol 
radiation  damping  varies  significantly,  depending  on  the  wave  speeds  assumed 
in  the  backfill  and  surrounding  media.  In  order  to  establish  a reasonable 
upper  bound  on  radiation  damping,  it  was  assumed  that  the  initial  tangent 
I'onstrained  modulus  of  the  backfill  was  two  times  the  value  recommended  by 
WKS  on  the  basis  of  1 .ibor.itory  tests.  This  produced  damping  due  to  radiation 
effects  o'  about  21  percent  of  critical;  evaluation  of  the  percentage  is  some- 
what subjective.  A friction  model  was  developed  which  added  about  5 percent 
to  the  radiation  damping,  leading  to  a total  of  about  2b  percent  of  critical. 
The  model  used  in  the  calculations  probably  overestimates  the  damping  due  to 
friction  effects  because  signs  and  magnitudes  of  friction  forces  .are  based  on 
absolute  rather  th.iu  relative  velocities.  Hysteresis  in  the  soil  was  not 
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icins  idcrinl  quant  i t at  i vi' I y , but  it  is  lu' 1 i cved  capablo  of  C('nt  ri  but  inj;  Ifss  tlian 
3 to  ')  piTcont  . Structural  damping  was  not  ronsidorod. 

An  attempt  was  made  to  evaluate  damping  from  the  physical  test.  This 
li'd  to  larger  values  of  damping  than  appear  reasonable  and  was  considered  unsuc- 
I ess  f u 1 . 

Unlike  the  case  of  the  arch,  forced  vibration  of  rectangular  structures 
produced  resonance.  This  is  partly  explained  by  the  back! ill  conditions,  which 
iliffered  significantly  between  the  two  types  of  struttures.  The  arch  had  a com- 
pacted backfill  of  native  soil  whose  uniaxi.il  stress-strain  curve  was  initially 
concave  to  the  strain  axis.  In  contrast,  the  rectangular  structures  had  sand 
backfill  which  was  placed  by  raining.  The  uniaxi.il  st  ress-st  r.i  i n curve  for 
the  sand  w.is  initi.illy  convex  to  the  strain  .axis,  which  le.ads  to  lower  seismic 
w.ive  speeds  than  for  the  .arch.  The  p.ar.ametric  series  of  calculations  on  the 
•arch  shows  that  its  radi.ation  damping  is  gre.iter  than  for  the  rect.angular  struc- 
tures .and  could  be  .about  two  times  gre.iter.  The  result  of  the  different  back- 
fill conditions  is  th.at  the  rectangul.ar  structures  were  closer  to  iji  y;iciuio 
conditions  than  was  the  arch. 

It  is  concluded  that  effective  damping  can  be  increased  from  10  to 
12  percent  critical  to  2f)  to  25  percent  critical  and  possibly  higher  by  increas- 
ing the  wave  speed  in  the  b.ackfill  .and  surrounding  soil  by  .i  f.actor  ot  three  to 
four.  Damping  due  to  friction  is  at  most  one-fifth  that  dui'  to  radiation  effects 
in  the  present  study.  In  cases  where  the  .active  earth  pressure  due  to  overburden 
is  less,  the  influence  of  friction  is  likely  to  be  less.  It  is  concluded  further 
that  .1  large  part  of  the  difference  between  the  responses  of  the  arch  .ind  the 
rect.angular  structures  is  due  to  the  different  backfill  materials  .iiid  methods 
of  placement. 
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RhCO.'tMIA'DAT  ION’S' 

llif  tolli'wiio;  iiTiimmiiulat  ioTis  .iri-  diviiltil  into  i .it  t^tor  i es  ot  .in.il- 
vsis.  .iml  clovi  lopmiiit  of  sinulo  iloKroo-ol  - f rcoilotn  modi' 1 s ol  compU-y 

riict  . 

1.  An.ilvsis  — I ho  nodo  1 of  iho  .inli  with  trillion  and  hi^h  w.ive 
spoi'd  h.irkfill  (caso  t,  Tablo  J-A)  should  bo  tisod  to  simiil  ato  portion  ol  t ho 
frequonov  swoop  in  whioh  qiias  i -st  oidv  stato  oonditions  ari'  roaolu-d  at  a niimbor 
•it  I 1 oso  1 v-spaood  t n-qiii  III' i os . This  woiihl  olarilv  whotin'r  rosonaiuo  is  appar- 
ont  wtion  radiation  and  friction  ot  foots  aro  oyplioitly  imludod  in  t ho  modi- 1 
rathor  than  boiriK  lumped  as  in  modal  damping-  Thoir  valuiss  should  In'  t hi'  m.ixi- 
mum  that  can  in-  ro.isonably  assum'd. 

2.  lostin^ — futuro  tosts  on  larp,o  model  'Uructuros  should  hi'  por- 
formod  with  t ho  samo  t horounhnoss , including  soil  tostinj;  in  t hi-  hacklill  and 
vibrations  at  difforont  li'vols  of  maximum  loroo,  as  wore  tho  KSSK.'<-V  to.ts. 

Iho  possibility  ol  moasurins  rolativi'  d i sp  1 acomont  s at  tin.'  so  1 1 -s  t ruo  t uri' 
intorfaco  tan^tontial  to  tin-  structuro  surtaco  should  hi'  i nvost  i ^a  t od ; this 
would  indicate  whothor  friction  should  ho  includi'ii  in  tho  calculations.  It 

• ippoars  that  much  inf orm.it  ion  miHht  ho  also  dorivod  I rom  sm.illor  scale  model 
ti'sts;  .1  series  ot  parallel  vibration  tosts  .md  slunk  lo.uiinp,  tosts  in  thi' 

Blast  l.o.id  ta'iior.itor  would  show  whothor  thi'  fori  oil  vibration  tests  include 
till'  samo  damping  phenomena  that  are  pri'Sont  in  shock  ii'S(ionse.  lin.illy, 

.It  tent  ion  should  be  >;ivon  to  a free  vibration  ph.isi'  ol  I lie  tests  in  which 
damping  could  bo  estimated  by  ohservinn  the  loj;  decrement  ol  the  .ini|>  I i tude . 

5.  Development  of  SDOK  Models--if  iiu  orporat  iny,  20  to  10  percent 
critic.'il  d.impiny  in  .SIX)!-'  models  h.is  a significant  impact  on  target  ing  decisions, 
new  methods  of  I'valuating  the  e<iuivalont  fraction  of  critii'.il  viscous  d. imping 
from  test  data  on  complex  structures  should  be  devi'loped.  Ihose  who  usi'  such 
models  should  become  .aware  that  backfill  conditions  can  afloct  damping  signifi- 
cantly and  that,  if  they  do  not  know  the  h.ickfill  londitions  in  the  structures 
they  are  analyzing,  they  should  assume  a range  ol  damping  values  from  about 
12  to  1^  percent  to  2')  to  ’K)  percent. 
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Al’IM  NUIX  I 

MODI'  SIIAPKS  or  ARCH  Sl'RUCTl'RK  WITH  AND  WITHorT  HACKKIll. 

N.'itur.'il  f r<.qui‘iK' ii'S  and  mode  sliapi'S  of  the  arcH  structure  were  cal- 
•ul.ited  with  .ird  without  h.ickfill  in  I'lie  soil  .ind  b.ickfill  properties 

Were  based  on  secant  moduli  up  to  10  psi.  Although  it  was  sul)set)uent  ly  found 
that  using  initial  tangent  moduli  improved  tt\e  dyn.iniic  simulation  using  direct 
i n t eg,  rat  ion  metliods,  tin  miula  I extraction  an.ilyses  were  not  repeated.  Previous 
experience  ( Kt  f . 9)  indicates  tliat  tile  frequencies  would  not  be  s i gn  i f i cant  1 y 
ini  reaseii  (less  than  10-15  percent)  by  basing  the  soil  stiffness  on  tangent 
moduli  .ind  that  the  mode  shapes  would  not  be  clianged  at  all.  An  end-on  eleva- 
tion view  of  the  slructute  and  surrounding  soil  is  sliown  in  figure  2-1.  for 
the  nod.al  analysis,  ali  subsurface  boundaries  ot  the  soil  isl.ind  were  fixed  ex- 
cept nodes  on  the  plane  of  symmetry  (par.alU'l  to  t lie  pi. me  shown  in  Figure  2-1), 
which  were  lonstr. lined  to  move  in  the  plane.  i'he  mode  shapes  of  tlie  structure 
are  stiown  with  soil  stripped  away.  Figure  1-1  shows  the  undefc'rmed  structure. 
Figures  1-2  tlirougli  I-ll  sliow  mode  shapes  ol  the  uncovered  structure.  Figures 
1-12  throegli  1-21  show  mode  shapi’S  of  the  covered  structure. 
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Al’PKNDIX  II 

Rv\I)lATl(lN  DAMPINC  ASSOC lATKD  WITH  DPI’ORMATIONAL  MODICS  OP  A 
CIRCPPAR  CYI.INDPR  PMHPDDKD  IN  AN  INPINITK  KI.ASTIC  DOMAIN 


It  is  shown  in  SecLion  3 Lliat  radiation  damping  associated  with 
vihration  of  the  covered  arch  increases  significantly  as  tlie  stiffness  of  the 
siirroiinding  soil  increases.  It  appears  tiiat,  in  tlie  range  of  uncertainty  witli 
wliicli  wave  speeds  in  tlie  surrounding  soil  are  known,  significant  variation  in 
r.idiation  damping  can  occur.  Tliis  finding  suggests  tliat  backfill  conditions, 
sucli  as  degrees  of  compaction  and  contact  with  the  structure,  may  need  to  be 
considered  in  developing  dam|)ing  parameters  for  tlie  idealized,  single  degree- 
of-freedom  analytic  models  used  in  targeting  analysis. 

To  confirm  the  trend  obtained  in  the  finite  element  calculations, 
Weidlinger  Associ.ites  requested  that  another  DNA  contractor,  l.ockheed  Paio 
Alto  Research  l.aboratory,  apply  a closed-form  analysis  technique  (Ref.  11) 
to  obtain  free  vibration,  mod  . 1 response  of  an  infinitely  long  cylinder  em- 
beddc'il  in  an  infinite  elastic  domain.  The  principal  investigator  at  l.ockheed. 
Dr.  T.  I..  (leers,  supplied  the  responses  for  tliree  problems  in  which  the  shell 
par.ameters  were  held  constant  while  the  wave  speeds  in  the  soil  varied  in  the 
ratio  1:2:A.  i'he  geometry  and  material  properties  of  the  cylindrical  shell 
and  soil,  di'fined  by  following  non-dimensional  ratios,  are  roughly  the  same  as 
those  of  the  arch. 


.0017 


= -^  = 1/2 

s c 

P 

whei  e 


h,  .1  “ thickn<‘S.s,  radius,  respectively,  of  the  shell 
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p , p = mass  density  of  tlie  shell  and  soil  medium,  respectively 
n m 


c , c = shear  and  dilatational  wave  speeds,  respectively,  in  the  soil 
s p 

medium 


The  wave  speeds  in  the  soil  are  defined  througli  a further  dimensionless  ratio 


o c 

V 


where 

K 

? o 

c = bar  velocity  in  the  shell,  c = 7,,  where 

o ■ o p ( 1 - V ; 


b , V are  Young's  Modulus  and  Poisson's  ratio,  respectively,  for 
the  shell  material. 


In  the  three  cases  considered,  B was  varied  by  varying  c the  dilatational 
wave  speed  in  the  soil.  The  range  of  parameters  was  as  follows: 

e P' 

Ca^e  _o 

1 47.6  (softest  soil) 

2 2'3.8 

1 11.9  (stiffest  soil) 

Tile  Lime  histories  of  radial  displacement  .ire  presented  in  Figures 

Il-l  through  11-4  for  modes  n = 0 through  n = ) ior  .an  impulsive  excitation  ol 

the  e.xtensionaj  modes.  The  r.idial  d i sp  1 ticemen t w for  mode  p and  the  time  t are 

normalized  to  a and  c /a,  respectively.  The  norm.i  1 iz.at  ion  parameter  for  time 
P 

makes  it  appe.ar  in  the  figures  .as  if  the  natural  f requeue  ii's  of  the  shell  are 
tirastically  reduced  if  the  shell  is  placed  in  stiffer  soil.  This  is  mainly  an 
illusion  c.aused  by  the  nond  imens  lon.a  1 p.ir.ameter,  as  can  be  seen  by  dividing  the 
time  sc.ale  in  cases  2 and  ) in  Figute  Il-l  by  2 and  4,  respectively.  Some  re- 
duction of  natural  frequency  occurs,  however,  app.arently  due  to  large  effective 
damping  ratios;  see  c.ase  "1,  Figure  II-4. 

'File  m.ain  purpose  of  performing  these  ca  1 cti  1 at  i ons  is  to  see  how  the 
equivalent  fraction  of  critic.al  viscous  damping  varies  with  wave  speed  in  the 
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soiJ.  In  these  cases,  damping  is  entirely  of  the  radiation  type.  There  is 
no  ambiguity  in  evaluating  the  log  decrement  for  eacli  displacement-time  his- 
tory because  each  represents  tlie  oscillation  in  one  mode  only,  and  the  decay 
rates  for  deformational  modes  are  negative  exponentials.  The  fractions  of 
critical  damping  are  sliown  in  Table  II-l.  Tlie  table  sliows  a trend  of  higher 
damping  with  increasingly  stiffer  soil,  which  is  the  same  trend  exhibited  by 
the  finite  element  calculations  of  the  arch.  When  the  properties  of  tiie 
system  are  fixed,  the  amount  of  damping  is  higlily  dependent  on  tlie  mode  number. 
For  example,  in  case  2,  mode  1 is  lightly  damped,  mode  2 is  about  critically 
damped  and  mode  4 is  liighly,  yet  suber i t ical 1 y , damped.  Tliis  clearly  illus- 
trates the  frequency  dependence  of  structure-medium  interaction.  The  n = 1 
mode  involves  rigid  body  motion  only  and  hence  the.  norma] ixed  responsL  appears 
to  be  independent  of  wave  speed  in  the  soil.  Tlie  damping  for  this  mode  is 
logarithmic  in  nature,  unlike  that  for  the  deformational  modes  wliich  is  expo- 
nential. Damping  in  the  rigid  body  mode  is  very  large  and  no  oscillation 
occurs . 
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Figure'  11-2.  Norm.il  i zfil  Modul  Dlspl.u'umeiU  for  ii  = 1 Modi*  Vor.sus  Normal  iz«'d 
Tinif  for  a Stiell  Fmboddod  in  Tliroo  Difforenl  Flnstic  Soil.s. 
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3.  Norm.i I i zed  Mod.il  Displareraent  for  ii  » 2 Mode  Versos  Normal  izeti 
Time  for  a Shell  I!mbedded  in  Three  Different  Klastic  Soils. 
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Fi^;iirf  I 1-4.  Normal  izeit  Modal  I)i  splacomenl  for  n = 3 Modo  Versus  Normalized 
Time  for  a Shell  Kmbeddcd  in  Three  Different  blast  ii-  .Soils. 


1 


J 


TABLE  n-1. 

EQLIVALEN'T  FtLACTIONS  OK  CRITICAL  VISCOUS  D/V-IPINC  FOR  INFINITELY  LONC 
CYLINDER  IN  INFINITE  SOIL  MEDIUM  SUBJECTED  TO  IMPULSIVE  LOADINC 


Mode 

n 

I 

(sof t soil) 

Case 

__2_ 

3 

(stiff  soil) 

1 

0 

.06 

. 1 3 

.24 

2* 

1 

>1.0 

>1.0 

>1.0 

i 

2 

.61 

a-KO 

M .0 

4 

3 

.19 

.41 

.51 

*Modi'  2 (n  “ I)  is  .1  rl^td  Body  t riins  I at  ion  mode  for  which  damping  is  large. 
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Al’l>i:\’DIX  III 


CHi'CK  I'ROlU.KMS  KOR  IMPLEMKNTATiON  (IF  FRICTION  IN 

THRKli-DIMFNSIONAL,  LINEAR  FINITE  ELEMENT  CODE 

riiL'  implomentat  ion  of  a friction  model  in  GENSAP  required  numerous 
coding  changes,  and  two  ciieck  problems  were  devised  to  tost  the  new  model. 

Tliese  are  illustrated  in  Figure  III-l.  Tlie  simulation  of  the  block  sliding 
on  the  plane  is  the  same  problem  for  which  Den  Hartog  (Ref.  2)  derived  equiva- 
lent viscous  and  Coulomb  damping.  The  friction  coefficients  in  the  expression 

f = pN 

was  chosen  such  that  p = 0.20.  This  produces  10  percent  equivalent  critical 
viscous  damping.  Initially,  a force  having  rise  and  decay  times  each  equal 
to  a quarter  of  the  natural  period  of  the  system  was  applii'd  to  the  block  to 
set  it  in  motion.  Disp lacement - t ime  histories  for  the  block  with  and  without 
Irictlon  are  shown  in  Figures  iIl-2  and  III-3.  Eliminating  the  damping  due  to 
the  integration  algorithm,  which  accounts  tor  the  decrease  in  amplitude  in 
Figure  III-2,  indicates  that  the  friction  is  responsible  for  approximately 
8 percent  d.amping. 

A more  ambitious  check  solution  was  then  obtained  for  the  two- 
dimensional  arch  subjected  to  a point  load  at  the  crown,  which  is  illustrated 
in  Figure  Ili-1.  The  di.iensions  of  this  arch  and  the  three-dimensional  one 
of  ultimate  interest  are  identical.  However,  the  mass  and  stillness  properties 
of  the  soil  were  I'mltted  from  the  two-dimensional  example.  In  this  calculation, 
the  frictional  forces  K were  computed  separately  tor  each  node. 

F = (pa^^)A 

where 

= the  normal  stress  at  each  node,  respectively,  due  to  overburden 
st  ress 

The  value  of  n varies  liue  to  variable  depth  of  a hypothetical  overburden  an  I 
n 

to  varying  angle  between  the  normal  to  the  arch  and  the  vertical. 
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Anolher  ilelnil  in  imi)lt'menL  ing  tiie  model  is  to  decide  on  tlie  direction 
ol  the  friction  force.  It  must  lie  along  the  tangent  to  the  arcti  at  each  in)de, 
liut  there  are  two  possible  directions.  Tlie  choice  is  made  by  comparing  the  direc- 
tion of  the  velocity  vector  with  the  normal  and  requiring  the  friction  force  vector 
to  lie  on  the  oppo.site  side  of  the  normal.  This  is  illustrated  in  Figure  IIT-1. 

As  in  the  other  example,  the  effect  of  damping  is  evaluated  by  apply- 
ing an  impulsive  force  to  the  crown  and  noticing  the  decrement  in  amplitude  of 
tlie  crt'wn  with  successive  cycle.s.  The  cases  shown  in  Figures  1II-4  and  111-5 
have  i.  = 0.4  and  zero  friction,  respectively.  As  in  the  case  of  the  block  slid- 
ing on  a frictionless  plane,  there  is  significant  decrease  in  amplitude.  In  this 
case,  it  is  partly  due  to  the  integration  algorithm  and  partially  duo  to  tlie 
nature  of  the  dynamic  response.  The  decri'ment  in  the  case  with  friction  is  much 
greater . 
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I>epartment  of  the  Interior 
bure.iu  of  Mines 

ATTN:  Tech.  lib. 

DfPART>fKN7  OF  DLUNSK  CONTRAt.loRS 

Aerospa.  e C.oi  poral  ion 

ATT.N:  le.  b.  Info.  Services 

ARbablan  Assoi  lates 

ATTN:  M.  AKbablan/F.  Safford 

Applied  Theory.  Inc. 

2 cy  ATTN.  John  t.  Triiiio 

Avi o Reaearrh  A Systems  Croup 

ATTN;  Research  Lib..  AH  JO,  Km.  7201 

Hattelle  Memorial  Institute 


ATTN; 

Tec hnlral 

Library 

Hie 

BDM  Corporat ion 

ATTN; 

Tec  hnlcal 

Library 

The 

BoeInK  Company 

ATTN; 

Acruspac  e 

Library 

I IT  Resear,  h Institute 

ATTN:  Milton  R.  Johnson 
ATTN:  R.  1 . We  Ich 
ATTN:  le>  iinic.il  Library 

Institute  for  Defense  Analyses 

ATTN:  IDA  Mbr-irian,  Ruth  S.  Smith 

Kaman  Avi Dyne 

Division  of  Kaman  Sciences  Corp. 

ATTN:  Tinhniral  Lib.  ary 

ATTN:  i.  S.  Cr  I s.  tone/N.  P.  Hobbs 

K.iman  Si  iem  «‘S  Corporal  ion 
ATTN:  Library 

Locktieed  Missiles  A Sp.it  <>  (Company.  Ini. 

ATTN:  Te.  h.  Inlor,  Center,  D/COLL 
ATTN:  Tom  (,eers,  D/S2-li,  Bldg.  20S 

laiveJai  e Fixin.l.ir  ion  for  Medl.al  hJu.ation  A Rsi  fi. 
ATTN:  Tethni.al  Library 

ATTN;  Asst.  Dir.  of  Res.,  Robert  K.  Jones 

Ml  Donnell  Doujil.is  Curpor.it  ion 
ATTN:  Robert  W.  H.ilprln 

Merritt  (sses.  Ini*orpi'r.»ted 

ATTN:  TechnUal  Library 

The  Mitre  Corpot.it  Ion 
ATTN:  Library 
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DFPARTMFNT  OF  DFFKNSK  (dNTRi\rTORS  ((.ont  inutnl) 


DFI'ARTMKNT  OF  DFFF.NSF  rONTRAC TORS  a.t.ni  I noiM  ) 


S.»th.in  M.  Nrwmark 
(x^nstiltln^  Khk  Inofr  ln>«  Servlt'cs 
BlOfeA  Civil  FtiK  IneiT  I ng  Hulldln>; 
I'nivrrMliy  ot  Iinn(.'l8 

ATTN:  Nathan  M.  NouTiwirk 

Physics  lntvrn.it  lonal  Company 


ATTN 

Doc  . 

Con . 

for 

F;.  T.  Mr«>re 

ATTN 

I)t>c . 

(on . 

for 

Dennis  Orph.il 

ATTN 

IVc  . 

Con . 

for 

l..irry  A.  Behrm.ann 

ATTN 

Doc . 

Con. 

tor 

Robert  Swift 

ATTN 

Iku  . 

Con. 

for 

Tech.  Lib. 

ATTS 

IkK- . 

Con. 

for 

(.harlcK  (.tulfrev 

R & f)  Ahk«)i' iat  t'S 

ATTN;  William  B.  WrlKhl,  Jr. 

ATTN:  Hvnry  Cooper 
ATTN:  Albert  L.  Uitter 
ATTN:  Harold  L.  Brodc 
ATTN:  J.  0.  Lewis 

ATTN:  Tvelinl.al  Library 
ATTN:  Reberi  Port 

Sciente  Appl  1«  .It  ions,  Inc. 

ATTN:  D.ivld  Bernstein 
ATTN;  D.  F.  Maxwell 

St  ieni  e Applications,  Inc. 

ATTN:  Technical  Llbr.iry 

Soul  liweHt  Re  Heart  h Inst  i t ut  t* 

ATTN:  A.  B.  Wenzel 

ATTN:  Wilt  red  E.  Baker 

Stantord  Researrh  Institute 

ATTN:  (ieor^e  R.  Abrahamyon/ B.  R.  (Listen 

Systems,  Science  It  Software,  liu  . 


ATTN: 

Thom.is  D.  Rini‘y 

ATTN: 

iKtn.ild  R.  (>rine 

ATTN; 

Ttuhnli.il  Library 

ATTN: 

Ted  Cherry 

Terr.i 

Ti'k. 

Iiu  . 

ATTN: 

Sidney  (ireen 

ATTN: 

1 ei  hf) itM  1 Library 

Tet  ra 

Tvi  fi. 

, 1 nc . 

ATTN; 

LI -San  Hwang 

ATTN: 

let  fin  it  a 1 Li  br.iry 

TRW  Defense  L Sp.ice  Sys.  Croup 

ATTN:  Pr.ivin  Hluilt.i,  RI-IIO'. 

ATTN;  Tech.  Into.  ( enl  er / ■ 1 'MO 
2 < y ATTN:  Peter  K.  Dal . HI /^1 70 

FRW  Defense  & Sp.u f Sys.  (.roup 
San  Bern.irdino  oper.ai  ions 

ATTN:  1.  Y.  Wony,  S27/71? 

Iniversal  An.ilyiics,  In«  . 

ATTN:  f.  1.  Field 

I'RS  Researt  h (T»mp.my 

ATTN:  TetlinIc.iI  I.ihrirv 

The  irlc  H.  W.my  ( tvll  l.tt^  i(«eer  inc.  Rs»  h.  I i-  , 
I’niversitv  Slat  it)n 

ATTN:  Neal  Haum 

Washlnp.lon  St.iie  Cnlversitv 
/Vim  In  I St  rat  Ive  Ory,.in  l/at  lop 

ATTN;  .\rthni  Miles  Moht'rf  ;or  'eorce  Duval 

WeltlliiiKer  Assoi  . (.onsultlnK  luKim-ers 
ATTN:  Melvin  L.  Bart*n 

Weitllinxer  Assot  . (airiHultlnR  Inylnetrs 
'«  cy  ATTN:  I.  Isenhery. 

AT  IN:  H.  S.  l.evine 
ATTN;  S.  H.  Pany 

West  itlKhoUKe  Fleet  ri(  Corp. 

Marine  Division 

ATTN:  W.  A.  Volz 


I 


